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O 8 Abstract 
<N 

(3J[) 9 The effect of a spatially inhomogeneous static magnetic field on the characteristics of capacitively 

<^ io coupled radio frequency discharges is investigated, with a focus on the sheath dynamics and the 

11 energy distribution functions (IEDFs) of the ions impinging the electrodes. Numerical studies 

12 are performed with two different kinetic algorithms. The first is a Particle-in-Cell (PIC) code 

13 which takes into account the entire discharge; the second is the novel Ensemble-in-Spacetime 

14 scheme (EST) which resolves only the sheath. When EST is fed with the sheath voltage and ion 

15 flux calculated by PIC, excellent agreement of the IEDFs is found: the applied magnetic field 
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16 predominantly influences the plasma bulk, i.e., the electron density and the ion flux to the walls; 

17 the dynamics of the sheath itself is only indirectly affected. It is concluded that EST may be used 
is as an efficient postprocessing tool to obtain the IEDFs in cases where only simplified (i.e., global 
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19 or fluiddynamic) information is available. 
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20 I. INTRODUCTION 



21 In the past three decades plasma processing such as surface activation, deposition, or 

22 etching has become of great importance for a large number of technological applications. 

23 With the goal of optimizing plasma processes, several types of plasma sources were been 

24 successfully introduced and studied, i.e., hybrid plasma sources, capacitively coupled ra- 

25 dio frequency discharges (CCRFDs) driven by more than one frequencies, or magnetized 

26 CCRFDs which are used for magnetically enhanced reactive ion etching (MERIE). It has 

27 been recently shown that dual frequency CCRFDs which use the electrical asymmetry effect 

28 (EAE, [TJ [2]) are most suitable to allow for a separate control of the ion flux and energy at 

29 substrates (or electrodes), which is one of the main goals in the context of plasma process- 

30 ing. In particular, the ion energy distribution function (IEDF) is one of the most important 

31 key parameters [3j 0] . By utilizing the EAE it is possible to invert the electrode function- 

32 ality, i.e., to change the DC self bias voltage from a large negative value to a large positive 

33 value even for geometrically symmetric CCRFDs with the driven electrode area equal to the 

34 grounded area. It has been further shown that geometrically asymmetric CCRFDs (small 

35 driven electrode and large grounded area, or vice versa), which produce a natural DC self 

36 bias voltage, can be completely symmetrized, in the sense that the natural DC self bias 

37 voltage can be completely compensated [3 [6] . 

38 In this work we study the effect of an externally imposed static magnetic field in com- 

39 bination with the EAE on the characteristics of dual frequency CCRFDs. We focus on the 

40 sheath dynamics and the IEDF at the electrodes. In contrast to the MERIE concept where 

41 usually a spatially homogeneous magnetic field is applied [7HTT]. we investigate the case of 

42 a spatially inhomogeneous magnetic field [12], which is able to (over-) compensate or even to 

43 amplify the EAE. 

44 The manuscript is organized as follows. After a short introduction of the discharge 

45 configuration and the operation regime, a short review of the employed numerical schemes, 

46 namely a one-dimensional Particle-in-Cell algorithm using Monte-Carlo collisions (PIC) and 

47 the Ensemble-in-Spacetime scheme (EST), is presented |13j . We discuss the simulation 

48 results for different magnetic field configurations focussing on the sheath dynamics and the 

49 IEDFs. We conclude that EST may be used as an efficient postprocessing tool to obtain the 
so IEDFs in cases where only simplified (i.e., global or fluiddynamic) information is available. 
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5i II. DISCHARGE SETUP 



52 In magnetized discharges (e.g., MERIE discharges) often a slowly rotating static and 

53 usually spatially homogeneous magnetic field is applied, with the goal to obtain a high 

54 electron density at very low gas pressures in the Pascal or even in the sub-Pascal range. The 

55 spatial dimensions of such discharges are typically of the order of a few centimeters, and 

56 the magnetic field strength can vary up to tens of mT. Under such conditions the electrons 

57 are magnetically confined while ions are not. The applied magnetic field predominantly 

58 influences the plasma bulk, i.e., the electron density and the ion flux to the walls, since the 

59 electrons suffer an increased number of collisions. It is important to note that the magnetic 
eo field does not affect the ion motion within the sheaths. Because of their large mass the 
ei Larmor radius for the ions is large compared with the typical length scales of the system. 

62 One can easily obtain a ratio of Tl and the sheath width s of r^/s = u P i/Qi ~ 100. Here 

63 ujpi and Q{ are the ion plasma frequency and the ion cyclotron frequency, respectively. 

64 In this work we investigate the discharge configuration depicted in figure [Tj In a geo- 

65 metrically symmetric plane electrode setup, with an electrode gap of L = 5 cm, an external 
ee magnetic field is applied near the left electrode. Within a small corridor (indicated by the 
67 dashed lines) the magnetic field can be considered parallel to the electrodes along the x- 
es direction. In contrast to MERIE setups where the magnetic field is spatially homogeneous, 
69 it is here assumed to decay from the left electrode into the plasma as 



70 The magnetic field amplitude is chosen to be Bo = 6 mT and the decay length is Ib = 5 

71 mm. The profile of B x (z) within the corridor is shown in figure [2j 

72 Further clarification is required to argue that the chosen one- dimensional approach in 

73 conjunction with the applied inhomogeneous magnetic field does not contradict Maxwell's 

74 equations. As mentioned above we use a one-dimensional approach in ^-direction to describe 

75 the discharge and apply a static magnetic field directed in the x-direction. If we assume 

76 that the magnetic field induced by the plasma current is negligible, Ampere's law given by 

77 Vx5 = /j, j requires 



B 



BqX 



(1) 




dB x dB z 



(2) 



dz dx 



3 



78 This necessarily leads to a physical magnetic field configuration as shown in figure [Tj This 

79 two-dimensional field configuration is similar to that of a "half magnetron" . Such a discharge 
so has of course to be treated by means of at least a two-dimensional simulation approach, 
si However, in this work we assume that the discharge can be described by a one-dimensional 

82 model. Such a simplified approach captures the essential physics because within the small 

83 corridor (indicated by the dashed line in figure [T]) the magnetic field lines are almost perfectly 

84 parallel to the electrodes. The steady state current in this configuration is therefore almost 

85 completely perpendicular to the electrodes so that the one-dimensional approach is justified, 
se In order to allow for the EAE a driving voltage V(t) of two consecutive harmonics, i.e., 
87 V(t) = 250 V(sin(2/T x 13.56 MHz t + Ay?) +sin(47r x 13.56 MHz t)), is connected through a 
as blocking capacitor of Cb = 1.5 nF to the left electrode. The relative phase A(f between the 

89 two harmonics is the control parameter for the EAE and varies between A<p = [— 7r/4, vt/4]. 

90 The right electrode is grounded. The pressure of the background Argon gas is fixed at p = 1 

91 Pa. 
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92 III. NUMERICAL ANALYSIS 



93 To analyze the magnetized (and non-magnetized) CCRFD described above, we employ 

94 two conceptually different kinetic approaches. The first is an explicit electrostatic ld3v PIC 

95 code (see e.g., [2]). Within a single PIC cycle (cf. figure [3]) the Poisson solver, the particle 

96 pusher and the Monte-Carlo collision module are successively iterated until convergence 

97 of the overall simulation is achieved. To couple the different modules, a first order field 

98 interpolation and charge assignment is performed. The particles are moved according to the 

99 Lorentz force using a push algorithm based on Boris' approach [T5]. Collisions of electrons 

100 and ions with the fixed background gas are performed in the frame of a slightly modified 

101 null-collisions method [16], with respect to the classical method proposed by Skullerud [T7] . 

102 Since Argon is used as the background gas the collision processes include elastic scattering 

103 of both electrons and ions, ionization and excitation due to electron collisions, and backward 

104 scattering (i.e., charge exchange collisions) of ions. Cross-section data are taken from Phelps 

105 et al. and the LXCat Database [T%] IT5]. 

we The second approach relies on EST which resolves only the sheath, unlike PIC which 

107 treats the entire discharge. EST is an iterative algorithm based on the solution of a set of 

108 kinetic equations for the ions, Boltzmann's relation for the electrons, and Poisson's equation 

109 for the electric field. Similar to PIC, EST employs the null-collision method to account for 
no elastic and charge exchange collisions of ions, but does not take into account ionization and 
in excitation processes. EST is fed by two input parameters: the first is the ion flux given 
n2 at the sheath edge. The second is the sheath voltage which can have an arbitrary periodic 
n3 waveform. Because of its efficient convergence behavior, EST allows for fast and kinetically 
114 self-consistent simulations of sheaths of CCRFDs. More details concerning the mathematical 
us description and the validation of the model can be found in [13J. 

lie In order to compare the results of the two different approaches, we extract the ion flux 

in and the sheath voltage from the PIC simulations. For this we calculate the sheath width 

us s(t) using the definition [20] 

ps(t) POO 

/ riidx = (th — n e ) dx. (3) 

JO Js(t) 

n9 Ion flux and sheath voltage are then used as input parameters for EST. 
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120 IV. RESULTS AND DISCUSSION 



A. Magnetized CCRFD 



122 We start our analysis by comparing two different discharge scenarios, without (case I) 

123 and with (case II) an applied static magnetic field at the left electrode by means of PIC sim- 

124 ulations. The simulation parameters are specified in section |TTJ The relative phase between 

125 the two driving voltages is set to Atp = 0. Figure [4] shows the space time dynamics of the 

126 electron densities for the two different cases. For case I we obtain, as expected, a symmetric 

127 density profile and a symmetric sheath dynamics. In contrast, in case II the discharge shows 

128 a significant asymmetry due to the applied magnetic field. One can observe that the mag- 

129 netized sheath in front of the left driven electrode is much smaller than the non-magnetized 

130 sheath in front of the right grounded electrode. This asymmetry can also be seen in figure 

131 [5] where the (time-averaged) ion density profiles are shown. It is clearly visible that the 

132 maximum is shifted and the slope is steeper towards the magnetized region in front of the 

133 left electrode. This asymmetry in the electro density (and the ion flux) is the reason the 

134 different sheath widths. Figure [6] compares the temporal behavior of the sheath widths for 

135 the two different cases of both the left (driven) and right (grounded) sheath. While in case 

136 I both sheaths show very similar oscillations, in case II the magnetized sheath is on average 

137 significantly smaller than the non-magnetized sheath (right). One can clearly observe a total 

138 collapse of the magnetized sheath at 10 ns. This magnetically induced asymmetry results 

139 in a substantial DC self bias voltage (of about 130 V). This can be seen in figure [7j where 

140 the time-averaged potential profiles for the two different cases are compared. 

141 Due to the high ion mass and a relatively low ion velocity, the magnetic field in front of the 

142 left electrode is not able to significantly confine the ions. The main effect introduced by the 

143 external magnetic field is the efficient electron heating near the left sheath. One can observe 

144 a locally increased plasma density and an increased ion flux to the electrode. Therefore, 

145 despite the magnetic field, we attribute the differences in the energy distribution functions 
we of ions impinging the electrodes (cf. figure [8]) mainly to the different plasma density profiles 
147 (thus the ion flux towards the electrodes) and to the sheath voltages accelerating the ions. 



148 The validity of this hypothesis is verified in section IV C, where IEDFs are calculated using 



149 the non-magnetic EST model. (Note that despite of the geometric and electric symmetry 



150 of case I, one may observe a small asymmetry in the sheath potentials and therefore the 

151 IEDFs. The reason for this is the difference of the sheath expansion and the resulting beams 

152 of energetic electrons caused by the corresponding electric fields in the sheaths. This effect 

153 has been previously reported by Schulze et al. [21]). 



154 B. Magnetized CCRFD with Electrical Asymmetry 

155 The parameter which controls the EAE is the relative phase between the two driving 

156 voltages Aip = 0. By introducing a spatially inhomogeneous magnetic field, which itself 

157 produces an asymmetry, one is able to amplify or (over-) compensate the EAE. As an exam- 

158 pie, we compare the previously discussed PIC results for case II with the results obtained 

159 for the same discharge conditions, but for relative phases of A(f = ±7r/4 instead of Aip = 0. 
wo Figure [9] shows the IEDFs for this scenario. Comparing the cases (cf. also figure [8]) one 
lei can observe that for Acp = — tt/4 the EAE can be used to amplify the magnetically induced 

162 asymmetry (in terms of the ion bombardment energy). Using a relative phase of Aip = ?t/4 

163 leads to an over-compensation of the magnetically induced asymmetry. 



164 C. EST: Magnetic Field Influence 

165 We have claimed above that the magnetic field does not significantly influence the dy- 

166 namics of the sheaths and thus plays only a minor role when calculating IEDFs. To justify 

167 this hypothesis we analyze the sheath dynamics using EST, which does not include any 

168 magnetic field effects. For this purpose we compare the IEDFs for the two cases without 

169 magnetic field (case I) and with magnetic field (case II) obtained by PIC with the results 



170 from EST. As described in section III we extract the ion flux and the sheath voltage from 



i7i the PIC simulation and use then as input parameters for EST. 



172 Figure 10 shows the space time dynamics of the electron density for both simulation 

173 approaches for case II (magnetized case for Aip = 0). The sheath dynamics is in very good 

174 qualitative agreement. Deviations can be observed particularly only near the sheath edge. 

175 This deviation is attributed to the different collision models used in PIC and EST: The 

176 latter does not account for ionization and excitation processes, while the PIC simulations 

177 allows for both ionization and excitation, and therefore efficient plasma generation. 



178 A second reason (which has to be proofed) is nonlinear electron resonance heating due 

179 to the excitation of higher harmonics in the plasma current and thus the plasma series 

180 resonance [2"2"H2"o'] . The higher harmonics are clearly visible in the sheath dynamics. Of 

181 course, this nonlinear sheath bulk interaction is not included in EST results, since EST 

182 resolves one sheath only. 



183 Figure 11 compares the IEDFs for the unmagnetized case I obtained using both PIC and 

184 EST. It is evident that both models are in excellent agreement. We can therefore conclude 

185 that it is well justified to use the "non-magnetic" EST for studying the effect of the magnetic 



186 field on IEDFs. Figure 12 shows the IEDF for the magnetized case II. The results are again 

187 obtained using both PIC and EST. Also for this case the results are in excellent agreement. 

188 Only minor deviations can be observed when comparing the magnitudes of the double peak 

189 structure. It is important to note that EST yields the same IEDFs as PIC, regardless of 

190 wether the discharge is magnetized or non-magnetized. The reason is, that both the sheath 

191 dynamics and the ion motion in the sheath are mainly driven by the electric field, and only 

192 weakly affected by the magnetic field. 



8 



193 V. CONCLUSIONS 



194 In this work we analyze the effect of a spatially inhomogeneous static magnetic field on 

195 the characteristics of capacitively coupled radio frequency discharges, with a focus on the 
we energy distribution function (IEDF) of the ions impinging the electrodes and the sheath 

197 dynamics. By means of a one-dimensional Particle-in-Cell algorithm using Monte-Carlo 

198 collisions (PIC) we show that geometrically symmetric capacitively coupled radio frequency 

199 discharges can be asymmetrized by applying a spatially inhomogeneous static magnetic 

200 field. This magnetically induced effect is similar to the electrical asymmetry effect (EAE). 

201 We further show that the EAE can be amplified or (over-) compensated. 

202 We find that the novel Ensemble-in-Space time scheme (EST) which resolves the sheath 

203 only, is able to provide IEDFs almost equal to the ones obtained from the fully self-consistent 

204 PIC simulation. Furthermore, we find an excellent agreement of the two models for both 

205 cases, the non-magnetized and the magnetized case. We conclude that EST may be used as 

206 an efficient postprocessing tool to obtain the IEDFs in cases where only simplified (i.e., not 

207 kinetic) information is available. Since EST allows for an arbitrary number of ion species, 

208 pressure, and composition of the background gas, it is a viable tool for the technology 

209 oriented computer aided design (TCAD) of realistic plasma processes [26] . 

210 Since the proposed discharge configuration is rather complex, a two-dimensional analysis 
2n of this problem is being performed and will be published elsewhere. Furthermore, the 

212 complex heating mechanisms of the magnetized CCRFDs due to the field-particle interaction 

213 are of great interest and require a further analysis which is beyond the scope of this paper. 
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FIG. 1. Schematic of the discharge configuration. A static magnetic field is assumed at the left 
electrode which falls off along the z axis according to ([!]) . There are two symmetric electrodes and 
the plasma is confined by a dielectric wall. The dashed lines indicate the region of interest for this 
work. 
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FIG. 2. Magnetic field profile inside the discharge. 



14 



Weight fields to 
particle positions 



s 



\ 



\ 



Calculate electric 
potential and field 




Weight particle 
charges to grid 



Calculate new particle 

positions and velocities 
1 

I 



MCC: Create or remove 
particles, adjust velocities 



FIG. 3. Iteration scheme of the Particle-in-Cell algorithm. 
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FIG. 4. Electron density obtained using the PIC method without (top, case I) and with (bottom, 
case II) an external magnetic field. 
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FIG. 5. Time-averaged ion density without (dashed, case I) and with (solid, case II) magnetic 
field. The maximum is shifted to the left towards the magnetized plasma sheath. 
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FIG. 6. Sheath modulation s(t) (non-magnetized and magnetized case, respectively) for the left 
(top) and the right plasma boundary sheath (bottom). 
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FIG. 7. Time-averaged potential along the discharge without (dashed, case I) and with (solid, case 
II) magnetic field. Due to the magnetic field a strong dc self-bias is established leading to a high 
ion bombardment energy (cf. also figure [8]). 
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FIG. 8. IEDFs for the two cases without (top, case I) and with (bottom, case II) an external 
magnetic field. With a magnetic field the discharge is strongly asymmetric. This is also detectable 
in the resulting IEDFs. 
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FIG. 9. IEDFs for the magnetized case for different relative phases of Deltaip. The magnetically 
induced asymmetry (Ay? = 0, middle) can be clearly ampfhfies (Ay? = — 7r/4, top) and also over- 
compensated (Ays = 7r/4, top)using the EAE. 
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FIG. 10. Electron density in the left sheath region obtained with the PIC method (top) and the 
EST model (bottom) for the magnetized case (case II). Despite the discrepancy in magnitude, 
there is good qualitative agreement regarding the sheath dynamics. 
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FIG. 11. IEDFs obtained for the non- magnetized case (case I) and A</? = using the PIC (top) and 
the EST model (bottom). There is excellent agreement between both models in the non-magnetized 
scenario. 
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FIG. 12. IEDFs for the magnetized case (case II) and A<p = 0. The PIC simulation (top) is 
performed with a constant inhomogeneous magnetic field, while the EST model (bottom) does not 
explicitly include the magnetic field, but is fed with input data obtained from the magnetic PIC 
simulations. 



24 



